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A light reflection technique is used to measure quanti- 
tatively the surface elevation of Faraday waves. The per- 
formed measurements cover a wide parameter range of driv- 
ing frequencies and sample viscosities. In the capillary wave 
regime the bifurcation diagrams exhibit a frequency indepen- 
dent scaling proportional to the wavelength. We also pro- 
vide numerical simulations of the full Navier-Stokes equations, 
which are in quantitative agreement up to supercritical drive 
amplitudes of e ~ 20%. The validity of an existing perturba- 
tion analysis is found to be limited to £ < 2.5%. 

PACS: 47.54. -hr 47.20.Ma 47.20. Lz 



I. INTRODUCTION 

A sound understanding of hydrodynamic pattern form- 
ing systems is based on a balanced interplay between 
experiment and theory, both analytical and numeri- 
cal. During the past this concept has led to consid- 
erable progress, especially in case of Rayleigh-Benard 
convection (RBC) or Taylor- Couette flow (TCP), where 
an amazing level of quantitative agreement has been 
achieved. Another famous example of pattern forming 
systems is the Paraday experiment: Surface waves on a 
free liquid-air or liquid-liquid interface are excited by a 
sinusoidal vibration of the fluid layer in vertical direc- 
tion. This system has the advantage of fairly short time 
scales in combination with a very rich bifurcation behav- 
ior. During the last two decades the focus in this field was 
on nonlinear pattern selection [^]-^ , secondary instabili- 
ties the transition towards chaos 0, droplet ejection 
lH and stokes drift ^. However, due to the paramet- 
ric drive mechanism the mathematical description of the 
Faraday experiment is more complicated as compared to 
RBC or TCP and the quantitative understanding is less 
advanced yet. For instance, it was not until recently that 
a rigorous linear stability theory had been developed, 
valid for viscous fluids and realistic boundary conditions 
I pof . Nowadays the predictions of the linear theory and 
the experimental results agree within a few percent p]| . 
On the non-linear level, however, the agreement between 
theory and experiment can at best be considered quali- 
tative: In the framework of a weakly non-linear pertur- 
bation analysis different primary surface wave patterns 
with quadratic, hexagonal etc. symmetry have been pre- 
dicted 1^ . Even though most of them have been found 
in experiments |^,||, empiric and predicted phase dia- 
gram reveal only a qualitatitve coincidence. Moreover, 



other experiments operate with a two- frequency drive sig- 
nal |l3|jl|], or at very shallow filling heights or with 
complex fluids |lq-|lq|. The resulting pattern dynamics 
becomes more complex and exotic structures like super- 
lattices or osciUons |l| appear. In some cases a 
qualitative understanding based on symmetry arguments 
could be obtained . 

With regard to this situation it comes as a surprise 
that a systematic quantitative investigation of the sys- 
tem's major order parameter, the surface elevation, has 
not been undertaken yet. Such data is indispensable for 
a veriflcation of any nonlinear theory. It is the aim of the 
present work to provide extensive experimental material 
in order to flU this gap. A measurement technique is pre- 
sented appropriate to quantify the surface elevation of 
Faraday ripples. Our measurements cover a broad part 
of the parameter space explored in recent experiments on 
pattern selection in the Faraday experiment. In case of 
low viscosity fluids our flndings are expected to compare 
with the perturbation analysis of Zhang and Vinals 
- at least at a weak supercritical drive. Furthermore, 
a closer comparison yields a reliable estimate of the va- 
lidity range of this approximation. Farther away from 
threshold of the instability, quantitative theoretical pre- 
dictions are not available yet. Here we provide a numer- 
ical simulation by means of a flnite difference scheme. 
This computation allows for 2-dimensional solutions of 
the full Navier-Stokes equations. They are used for a 
comparison with our measurements on line patterns. 



II. THE SYSTEM 

We consider a fluid layer of thickness h with a free sur- 
face vibrated in vertical direction at a drive frequency 
fl. In the frame of reference co-moving with the con- 
tainer the liquid is subject to a modulated gravity ac- 
celeration g + a sin {fit) , where g is the gravitational ac- 
celeration and a the amplitude of the drive. The fluid, 
considered as being incompressible, is characterized by 
the kinematic viscosity p, the density p, and the sur- 
face tension a. The hydrodynamic problem is governed 
by the Navier-Stokes equation, in which the modulated 
gravity enters as a parametric drive. The boundary con- 
ditions are free slip at the fluid-air interface and no slip 
along the walls and the bottom of the container. If a 
exceeds a critical threshold ad^^, h,v, p, a) the surface, 
being plane at the beginning, undergoes the Faraday in- 
stability and standing waves with a wavenumber k ap- 
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pear. Usually these waves organize themselves in form of 
regular patterns of different possible symmetries (lines, 
squares, hexagons ...). In containers with a large lateral 
aspect ratio (the ratio between container dimension to 
wavelength of the pattern) the pattern selection is geom- 
etry independent and solely governed by the nonlinear- 
ities in the equations of motion. In contrast, by chang- 
ing the lateral container extension one can manipulate 
the selection process. For instance, by reducing one side- 
length of a beforehand quadratic container, a line pattern 
(oriented parallel to the shorter side) can eventually be 
enforced in a situation where squares would prevail oth- 
erwise. In the present paper extensive use will be made 
of this geometrical selection feature. 

Generically, the time dependence of the Faraday mode 
bifurcating out of the undisturbed plane surface is sub- 
harmonic, i.e. the surface oscillates at half the drive fre- 
quency uj — Q/2. The standing wave surface profile 77(r) 
at onset of the instability can be written in the form 



1 ^ 

i=l 



+ CXJ 



iki 



(1) 



Here r = {x, y) abbreviates the horizontal coordinates. 
The set of Fourier coefficients {Cn} determines the sub- 
harmonic time dependence. Exactly at the onset of the 
instability, a = Oc, one has Cn — for even n, while 
the odd coefficients are the components of the eigenvec- 
tor related to the linear stability problem. The spatial 
modes are characterized by the wave vectors , each car- 
rying an individual amplitude Ai. These quantities are 
determined by the nonlinearities of the problem and - if 
appropriate - also by the container geometry. In princi- 
ple the ki can have any length and orientation but usu- 
ally they are supposed to be equally spaced on the circle 
|ki| = k. Then the number N of participating modes 
determines the degree of rotational symmetry of the pat- 
tern: A*" = 1 corresponds to lines, = 2 to squares, 
A'' = 3 to hexagons or triangles, etc. Patterns with > 3 
are no longer translational symmetric, we refer to them 
as quasi-periodic. 



III. PREWORK BY OTHER AUTHORS 

A. Experimental 

There exist a considerable number of measuring tech- 
niques appropriate to investigate the surface wave dy- 
namics. They range from contactin g p ermittivity mea- 
surements poll over optical systems [ pi[ including inter- 
ferometry, radar back scattering up to x-ray absorp- 
tion techniques Most of these procedures apply to 
gravity waves, i.e. long wavelength surface waves. Each 
of them has its own limitations. Either the amplitude 
is recorded locally in space, or the detector is too slow 
to resolve the temporal spectrum. One of the strongest 



restrictions for optical reflection techniques comes from 
the poor reflectivity of most fluid- air interfaces. Only a 
few percent of the incident light is reflected and an even 
smaller contribution is due to diffusive back-scattering. 
Since the latter is crucial for standard visualization tech- 
niques such as holographic photography or triangulation 
p^ , these methods often do not work properly. One 
might think of spraying the fluid surface with tracer par- 
ticles like club moss seeds but this in turn affects the 
surface tension or the surface viscosity in an uncontrol- 
lable manner. In a pair of very recent publications ^5 26 
the authors propose sophisticated methods based on a 
colored illumination or on an array of microlenses. A 
standard visualization technique used in several Faraday 
experiments is the light shadowgraphy [^PJ2^. A beam 
of parallel light with a diameter comparable to the con- 
tainer size passes through the fluid layer. The deformed 
standing wave surface pattern with its heaps and hollows 
acting like an array of collecting and diffusing lenses is 
mapped onto a screen. However, this technique is re- 
stricted to very small surface deformations for which the 
profile can be recovered by ray tracing. Otherwise caus- 
tics occur, which make this simple and effective method 
break down. 

An alternative technique introduced by Wright, Bu- 
dakin and Putterman ||28|] is based on intensity losses due 
to diffusive scattering. These authors use Polystyrene 
colloids to provide light scatterers within the fluid. The 
transversing light intensity, weakly modulated by the lo- 
cal layer thickness, is recorded by a high-sensitivity CCD- 
camera to reconstruct the surface profile. As an example, 
a snapshot is presented taken at a rather strong drive am- 
plitude within the turbulent regime. 

To our knowledge, previous quantitative investigations 
on near-onset Faraday waves are restricted to the low- 
frequency regime, i.e. below fl ~ 30Hz. Douady p9| ] 
presents an investigation, where the deflection of a thin 
laser beam directed onto a wave node is recorded. To 
prevent the nodes from a gradual drift, Douady uses the 
" rimfuU technique" to fix the position of the pattern rela- 
tive to the container boundaries. That way the structure 
experiences a mode discretization dictated by the con- 
tainer geometry. Data for the elevation maxima and the 
linear relaxation time r are provided for a silicone oil at 
a viscosity of lOcS and drive frequencies between 20 and 
30Hz. 

Several other authors report amplitude measurements 
in small aspect ratio experiments at driving frequencies 
between 5 and lOHz p0|-p2t. In this situation a con- 
siderable mismatch between frequency and wavelength 
may occur. The resulting detuning renders the primary 
bifurcation hysteretic. An additional problem in these 
experiments is the moving contact line between the fluid 
and the container. It may also generate a hysteresis or 
even induce an irregular time dependent wave dynamics. 
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B. Theoretical 



1. Numerical 



A couple of numerical simulations on Faraday waves 
is worth mentioning in the present context. Zhang and 
Vinals reduce the original full 3-dimensional hydro- 
dynamic problem to a set of 2-dimensional non-local 
equations formulated in terms of the lateral coordinates 
only. This is the outcome of the so called " quasi-potential 
approximation" , which considers the bulk flow as being 
potential (inviscid). Vortical flow contributions in the 
viscous boundary layer beneath the surface are accounted 
for by an effective boundary condition. Clearly this anal- 
ysis is only valid for deep layers {h oo) and restricted 
to the low dissipation limit i/k^/il <^ 1. The derived 
equations for the surface elevation •q{r,t) are numerically 
integrated by a pseudo-spectral algorithm. The principal 
concern of these simulations is to obtain the most pre- 
ferred pattern and to investigate its resistance against 
secondary instabilities. Stationary patterns with differ- 
ent rotational symmetries have been observed. In par- 
ticular, drive frequencies at the transition between grav- 
ity and capillary waves, r2/27r ~ 30Hz, give rise to the 
most interesting structures: Quasi-periodic pattern with 
a 10-fold rotational symmetry [Q. In the capillary regime 
squares or lines are observed depending on the viscosity 
of the sample fluid. 

A somewhat different numerical procedure has been 
used by Schultz et al. Q and Wright et al. H]. Both 
investigations are based on the Euler equations for invis- 
cid fluids. Dissipation is accounted for by a phenomeno- 
logical linear damping term introduced afterwards. The 
numerical procedures used are respectively a boundary- 
integral method and a vortex-sheet technique. That way 
the profile of very steep 2-dimensional Faraday waves 
is investigated. Their findings comprise dimpled crests, 
the formation of a plumelike shape, or the beginning of 
droplet ejection. 

Obviously, all of the aforementioned numerical work is 
restricted to the low dissipation limit and to large layer 
thicknesses. Indeed, we are not aware of systematic nu- 
merical simulations for viscous fluids on the basis of the 
full hydrodynamic equations and boundary conditions. 



2. An 



The first step towards a theoretical understanding of 
pattern forming systems is a linear stability analysis. 
This gives access to the threshold amplitude Oc, the crit- 
ical wave number fee, and the most unstable mode. For 
mathematical convenience it is popular to assume later- 
ally infinite geometries. Although it is known for a long 
time that the stability problem for Faraday waves can be 
approximately mapped to a parametrically driven pen- 
dulum [p6| a rigorous stability analysis for viscous fluids 



dates back until recently ||T^. Quantities evaluated by 
this method will be referred to in the following as "re- 
sults of the exact stability analysis" . 

A nonlinear analysis suitable to predict the selected 
surface pattern just above the primaryinstability has 
been presented by Zhang and Vinals |12). They start 
form their reduced 2-dime nsional set of quasi-potential 
equations (see Sec. [II Bl) and perform a perturbation 
analysis for small supercritical drive e = {a — Uc) / Uc ^ 1. 
The analysis can be understood as a double expan- 
sion in e and in the dimensionless damping parameter 
7 = vk^ /fl ^ 1. Since /i — > oo is assumed, vortical flow 
contributions coming from the viscous boundary layer 
along the bottom of the container are ignored. The pat- 
tern wavelength 27r/fc is approximated by the inviscid 
dispersion relation 



= gfc 



(2) 



In our experiments the damping parameter is 0.1 < 7 < 
0.2 for the low viscosity sample SI and lies between 1.3 
and 1.7 for probe S2. Accordingly, the wave number k 
evaluated by Eq. (||) for SI agrees within 1% with the 
result of the exact stability analysis but it is off by 20% 
for S2. In the framework of the weakly non-linear analysis 
the surface profile is represented by Eq. (|l|) with \ki\ = fee 
and Cn = for " even. The mode amplitudes Ai are 
governed by the set of the amplitude equations 



N 



(3) 



where r is the time constant of linear damping, also an 
outcome of the linear analysis. r(0) is the non-linear 
coupling function, whose dependence on the sample spec- 
ifications and the drive frequency is analytically known. 
Evaluating P at the angle increments 6ij between two in- 
teracting modes and yields the set of cubic coupling 
coefficients, which governs the pattern selection process. 
By computing the stationary solutions of Eq.^ in the form 
\Ai\ ~ TZn with i = l...n and 



\A,\^n 



N 



(4) 



at different symmetry indices A'^ one obtains the satura- 
tion amplitude TZi for line patterns, 7^2 for squares, etc. 
As outlined in Refs. ||l^;0 that pattern with the smallest 
free energy 



N 



J2 m,)\A^\A,\ 



(5) 



is to be selected. Broadly speaking, in the capillary 
regime investigated here, it is either squares (at low sam- 
ple viscosities) or lines (at higher viscosities). Wave pat- 
terns with a higher order rotational symmetry are not 
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found since they routinely posses a larger free energy. 
Since our optical reflection technique works best with a 
l-dimensional surface modulation, we make use of the 
above mentioned geometrical selection feature and per- 
form all experiments in containers of rectangular cross 
section. That way we enforce line patterns at any inves- 
tigated sample viscosity and under all operating condi- 
tions. Under the assumption that the saturation ampli- 
tude for lines is not appreciably influenced by the con- 
tainer geometry we are able to compare our data with 
the Hne solution (iV = 1) of Eq. pi 



r]{x,t) = Acos{kcx) sm{—t). 
Here the saturation amplitude is given by 



(6) 



A = (7) 
where we have abbreviated Ai by A and r(0°) by F. 

IV. EXPERIMENTAL SETUP 

The surface wave excitation is accomplished by a vi- 
bration exciter with a maximum force peak of 4670N. 
Details are described elsewhere . The acceleration sig- 
nal is computer controlled, its amplitude and harmonic- 
ity are stabilized such that fluctuations are smaller than 
±0.2%. We use a container with a rectangular cross sec- 
tion (length 150mm, width 50mm) in the form of a sta- 
dium (see Fig. |]). This side length ratio is sufficient 
to enforce Faraday patterns in form of lines for all in- 
vestigated sample fluids and under all operating condi- 
tions. To minimize disturbances originating from menis- 
cus waves, the form of the rim is designed as a "soft 
beach" , where the depth increases up to its maximum 
of 5mm on a length of 9mm giving an inclination angle 
of ~ 34°. The curved sides of the stadium also help to 
suppress meniscus waves due to destructive interference. 
The vessel was covered by a glass plate to avoid pollution, 
evaporation and temperature fluctuations. The temper- 
ature of the container (typically T = 30° C) is regulated 
with an accuracy of ±0.1°. Two different Dow Corning 
200 silicone oils are used as sample fluids. For a com- 
plete specification see Tab. L The choice of the filling 
height h = 3mni is large enough to guarantee that the 
finite depth correction to the dispersion relation (||) can 
be ignored, i.e. tanh(fcft,) > 0.995 for all measurements. 

The knowledge of the wave number k is crucial for 
the interpretation of the elevation amplitude (see below). 
Therefore, before each run of amplitude measurements 
we evaluate the wave number k of the pattern by pho- 
tographing the free surface illuminated by a diffuse light 
source (see Fig. 

To record surface wave amplitudes we use a laser beam 
directed vertically onto the fiuid surface (see Fig. |^) . The 



cross section of the beam is widened to a diameter be- 
tween 0.75 and 1.25 of a pattern wavelength. The light 
beam reflected at the standing wave surface pattern hits 
a diffusive screen mounted above the liquid-air interface. 
The shape of the reflected light pattern depends on the 
surface wave structure. In case of lines, a bright light 
streak with sharp edges occurs on the screen. Its length, 
oscillating with the frequency of the external drive, is 
recorded by a CCD camera and digitized. The largest 
length during an oscillation cycle yields the maximum 
surface slope dr/{x)/dx\max- According to the geometry 
shown in Fig. S one obtains 



dx^ 



tan(/3) 



1 s 
tan ( — arctan — 
^2 



(8) 



Deflection angles up to /3max = 40° have been exploited. 
Minor effects due to the refraction of the light beam 
by the glass cover are corrected for. As seen in Fig. ^ 
the light rays marking the tips of the streak may orig- 
inate from two neighboring elevation nodes. This and 
other errors in s together with the inaccuracies in k and 
d sum up to a relative systematic error in drj/dxlma-^ 
which does not exceed 10% for dive amplitudes e > 0.5% 
(see Fig. ^). In case the surface proflle at the moment 
of maximum elevation is approximately sinusoidal with 
wave number k the amplitude jyjmax can be deduced from 
(I) via 



V\ 



k dx 



(9) 



The light pattern on the diffusive screen is recorded by 
a CCD camera situated vertically above it. The pictures 
are evaluated by a home made processing software. To 
that end each image is binarized and the longest distance 
between any two points of the light streak is extracted. 
An automatic adaption of the binarization threshold is 
implemented to compensate for the decreasing light in- 
tensity of the streak during an amplitude ramp. 

The measurement of the bifur- 

cation diagram drj / dxlme^Ki^) runs as follows: Starting 
at a drive amplitude of £ < 0.2%, the computer performs 
an automatic ramp in steps of Ae = 0.065% up to the 
maximum Smax- For sample SI it is Emax = 20% but only 
Emax = 2% for the more viscous probe S2. The fluctua- 
tions of the drive acceleration vary from ±0.05% at low 
drive (a ~ Ig) up to ±0.2% at a ~ 20g. Between each 
increment the ramp is suspended for a waiting period of 
1 minute to allow for the system to equilibrate to the new 
situation. Then a series of 25 snapshots of the light streak 
is taken at regular intervals of 20s. This yields the av- 
erage streak length and the statistical error as indicated 
by the error bars in Fig. IJb. After each upwards ampli- 
tude ramp a second scan in opposite direction down to 
£ — —2% is performed to check for a possible hysteresis 
in the bifurcation diagram. 
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V. NUMERICAL SIMULATIONS 

In this paper we present extracts from our numerical 
simulations of the full nonlinear hydrodynamic problem 
adapted to treat 2-dimensional Faraday patterns in form 
of lines. A sketch of the implemented algorithm is given 
here, details will be presented elsewhere [Q. 

The line patterns are considered 2-dimensional in the 
x-z-plane, the y-dircction is ignored. The applied al- 
gorithm is based on a popular marker-and-cell-mcthod 
(MAC) 1^,0 extensively used to simulate thermal con- 
vection in fluids Therein the non-dimensionalized 
incompressible evolution equations read as 



condition at the bottom. For the present simulations pe- 
riodic boundary conditions in lateral direction are used, 
even though our algorithm allows to switch easily to a 
realistic no-slip situation. All simulations are performed 
by a mesh composed of 80 cells per wavelength 27r/fc at 
a time step size of 0.001. 

We emphasize that the present integration method 
does not suffer from the limitations of earlier algorithms 
(see Sec. IIIB 1 above). In particular, since it is based on 
the full Navier-Stokes equations, there is no restriction to 
the weakly dissipative limit. Moreover, our algorithm al- 
lows to study finite filling levels, even if the depth of the 
viscous boundary layer compares to the layer thickness 
h. 



dtu = -dxu^ - dziwu) + ^^jp- {dxxu + dzzw) - dxP , (10) 
dtw = -dx{uw) - dzw'^ + {dxxU + d^^w) - d^p, (11) 



dxu + d^w = (12) 



with u and w being the x- and z-component of the ve- 
locity, respectively and p being the pressure. To non- 
dimensionalize we have used the following unities: wave 
number k for length, $7/2 for time and pri^/(4fc^) for pres- 
sure. The above set of equations is solved on a staggered 
grid for u, w, p. The time integration is carried out by 
a forward time step, while the diffusive space derivatives 
and the pressure gradient is evaluated by central differ- 
ences. The convective terms are treated by a Donor- Cell- 
scheme. With a damped Jakobi-method the pressure is 
determined such that the incompressibility condition ( p^ ) 
is met. Since we are not interested in surface profiles 
with breaking waves or droplet ejection, r](x,t) is a sin- 
gle valued function. Its dynamics is determined by the 
kinematic boundary condition 



The dynamical boundary conditions 



P\z=n-'"--c' \z=n 



n= 4^ (1 + asin(2t)) 



f72 



(13) 



(14) 



(15) 



ensure the continuity of tangential stresses across the in- 
terface and the discontinuity of normal stresses due to 

,2 

the finite surface tension. Here a'^j = 2^^{V iUjV jUi) 
denotes the dimensionless viscous stress tensor, n{x, t) 
is is the surface normal vector in outward direction, and 
t{x,t) the tangential vector perpendicular to it. Note 
that the tangential condition is completely ignored in 
some previous free surface algorithms Here the im- 
plementation of Eq. ( p^ is accomplished by approximat- 
ing the discretized interface line by either horizontal, ver- 
tical or diagonal segments as suggested by Grieb Q . Be- 
sides the surface boundary conditions we impose a no-slip 



VI. RESULTS 
A. Measurements at low viscosity 

By virtue of the rectangular container shape the pri- 
mary Faraday pattern consists of lines oriented parallel 
to the shorter sidewall. This is in contrast to control ex- 
periments carried out in large aspect ratio vessels, where 
squares are the selected planform under otherwise iden- 
tical conditions. At the onset of the instability defining 
£ = the lines occur first in local regions of the surface. 
By increasing the drive up to e « 0.5% the line pattern 
spreads out over the whole surface. This non-ideal on- 
set is due to the spatial inhomogeneity of the drive. To 
rule out whether the finite longitudinal container dimen- 
sion gives rise to a mode discretization we scanned the 
drive frequency between ^1/2tt — 80Hz up to lOOHz in 
steps of 0.5Hz. Thereby the number of waves fitting into 
the container increased from 32 to 37 in a continuous 
manner. By virtue of the "beach like" container rim no 
stepwise behavior of the curves kc{^l) or ac{fl) could be 
detected. Moreover, the measured onset amplitude and 
wave number always agreed within 0.5% with the theo- 
retical prediction of the exact linear analysis computed 
for a laterally infinite fluid layer. 

A set of amplitude measurements performed at four 
different drive frequencies is shown in Fig. |^. Only 
the data obtained by up-ramping the drive are plot- 
ted since the corresponding down-ramps did not deviate 
significantly. There was no indication for a hysteretic 
transition. The highest drive amplitude achieved was 
£max — 40% (Fig. ||b) where the maximum surface eleva- 
tion reaches about 0.5mm, i.e. already 16% of the layer 
thickness. However in most runs the amplitude scan had 
to be stopped at £max — 20% due to an incipient defect 
dynamics. Over the whole investigated drive amplitude 
range, < £ < 20% (40%) no deviation of the non-linear 
wave number k{e) from its onset value fc^ could be de- 
tected. 

The bifurcation diagrams ?7|max(£) as shown in Fig. ^ 
demonstrate that a square-root like increase according 
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to the theoretical prediction 77 1 max = \/e/r (sohd hne in 
Fig. H) is limited to rather small drive amplitudes below 
£ ~ 2.5%. For a closer quantification we used the exper- 
imental data at < e < 2.5% and fitted the coefficient 
r. Good agreement with the prediction of Zhang and 
Vinals (see squares and solid line in Fig.^) is observed. 
It is interesting and to our knowledge unmentioned yet 
that the analytical expression for F/fc^ becomes ft- 
independent at large drive frequencies according to 



240 >^ 



(16) 



It can be seen from Fig. ^ that this asymptotics applies 
fairly well already at drive frequencies beyond f2/27r ~ 
60Hz. As a consequence of Eq. (^6|) the surface slope 
977/92: linax should bc asymptotically independent of the 
drive frequency. This is tested by Fig. ^, where the sur- 
face slopes associated with different drive frequencies ap- 
proximately collapse on a common master curve. This 
scaling persists even for e > 2.5%, i.e. beyond the valid- 
ity range of the perturbation analysis. 

Let us briefly discuss possible sources for discrepan- 
cies between pertubation analysis and experiment. As 
mentioned above, the appearance of line patterns aligned 
parallel to the shorter container side (transverse mode) 
is enforced by the rectangular container geometry. The 
result of this finite geometry effect may be twofold, linear 
and non-linear. The linear one is negligible as can be seen 
from the following argument: The container sidewalls 
provide a damping offset to both, the longitudinal and 
the transverse line mode. However, since the distance 
between the longer sidewall pair is smaller the threshold 
shift for the longitudinal mode is enhanced. Indeed, we 
even find the threshold of the transverse mode almost 
unaffected: The empiric Faraday onset Oc agrees within 
0.5% with the prediction of the exact stability analysis 
computed for a laterally infinite container. Clearly, this 
linear reasoning does not imply that the saturated nonlin- 
ear pattern amplitude remains unchanged too. However, 
estimating the geometry effect upon the coefficient F is 
a difficult task: To that end one had to redo the per- 
turbation analysis in terms of the container eigenmodes 
instead of the simpler plane waves. 

A second cause for deviations to the experimental re- 
sults is the restriction of the perturbation analysis to 
a pure sinusoidal time depeiidence as given by Eq. (|^). 
However, as outlined in Sec. HI B 2 the actual frequency 
spectrum of the Faraday mode is not monochromatic. In 
the investigated region the error in Oc resulting from this 
approximation lies between 10% and 20$. 

For a quantitative comparison with the experimental 
data at elevated drive amplitudes w e refer to the numer- 
ical simulations as described in Sec. [II B 1 and indicated 
in Fig. ^ by the open symbols. Up to elevated drive 
amplitudes of £niax — 20% we observe good agreement 
between simulation and experiment. The deviation is 
nowhere worse than 15% but considerably better at the 



higher drive amplitudes, where the systematic error of 
the detection method is smallest (see Figs. ^,d). Beyond 
predicting the local elevation maximum the simulations 
also provide access to the spatial anharmonicity of the 
surface elevation. We find that the spatial harmonics 2fc, 
3fc ... contribute to the Fourier spectrum by less than 
4%. This justifies a posteriori to equate the elevation 
amplitude ?7|,„ax with the ratio dr]/dx\max/kc as used to 
produce Fig. 0. 



B. Measurements at high viscosity 

Following the same procedure as in the previous section 
we performed a set of amplitude measurements on the 
more viscous fluid sample S2 (see Fig. ^) at a temperature 
of 30° C. Again the selected pattern consists of parallel 
transverse lines. However, unlike SI, this is the preferred 
planform in large aspect ratio containers too, being a 
result of the elevated viscosity of probe S2. Consequently 
the rectangular container geometry just determines the 
orientation of the lines rather than altering the selected 
pattern. 

Also in contrast to SI the more viscous probe S2 ex- 
hibits a wave amplitude, which grows rapidly with the 
driving acceleration (Fig. ^). The maximum deflection 
angle of /3max — 40° is reached at £ < 1%. Clearly, 
at those small drive all bifurcation diagrams show fairly 
well a square- root like increase according to Eq. (|^). 
However, the coefficient F as compiled from the data 
at < £ < 0.8% is an order of magnitude larger than 
for SI (circles in Fig. ^). Moreover, F deviates substan- 
tially from the prediction of the perturbation analysis 
(dashed line in Fig. But the latter observation is not 
surprising since both the low damping approximation as 
well as the infinite depth assumption are violated for S2: 
Note that 1.3 < 7 < 1.7 and the depth of the viscous 
boundary layer is 0.5mm, i.e. 20% of the layer thickness. 
Although the perturbation analysis does not quantita- 
tively apply here, we recover the same i7-independent 
scaling of the bifurcation diagram. The data for the 
slope Or] / dx\jna.x{£) collapse again on a master curve (see 
Fig. II) the same way as they did for SI in Fig. ||. 



C. Viscosity dependent measurements 

A final set of measurements is devoted to the viscos- 
ity dependence of the surface elevation. Fig. ^ shows a 
set of bifurcation diagrams obtained with S2 at different 
viscosities (see Tab. VH). This is accomplished by vary- 
ing the temperature of the probe. These measurements 
are performed at the drive frequency Q,/2Tr — 80Hz. In 
agreement with our previous observations the surface el- 
evation steeply rises as the viscosity increases. Fig. ^ 
shows the dependence of the coefficient F on the viscos- 
ity as fitted from the experimental data. By comparison 
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with the result of the weakly nonlinear analysis (solid line 
in Fig. ^) we conclude that the small damping approxi- 
mation 7^1 holds at best up to viscosities of ~ 50cS. 

VII. CONCLUSIONS 

We have presented a series of systematic amplitude 
measurements for stationary Faraday surface waves. The 
investigation is accomplished by a laser beam reflected at 
the oscillating surface. To facilitate the interpretation of 
the data the measurements are performed on line pat- 
terns, which are enforced by the rectangular container 
geometry. Due to the soft "beach like" boundary condi- 
tions a mode discretization is avoided. Bifurcation dia- 
grams of the maximum surface deflection vs. the drive 
amplitude are systematically recorded over a wide pa- 
rameter range of drive frequency and sample viscosity. 

The experimental data reveal that the perturbation 
analysis of Zhang and Vinals jl^ applies quantitatively 
to fluids with a viscosity of less than ~ 50cS and to very 
small drive amplitudes of not more than 2.5%. More- 
over, we observe that the surface slope scales almost in- 
dependently of the drive frequency. This finding is also 
supported by the analytical expression for the nonlinear 
coupling coefficient T as derived in the framework of the 
perturbation theory. Qualitatively this scaling behavior 
persists even up to drive amplitude of Emax — 20%, i.e. 
at operating conditions, where the perturbation theory 
is no longer applicable. 

For a quantitative comparison of our data at elevated 
drive amplitudes we provide a numerical simulation of 
Faraday waves on the basis of the full Navier-Stokes equa- 
tion. This new algorithm does not suffer from the stan- 
dard restrictions of the low dissipation limit and large 
filling thicknesses as used by other previous simulations. 
Good quantitative agreement with the empiric data is 
found up to the highest investigated drive amplitudes of 

Emax — 20%. 
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FIG. 1. Photograph of the container and a surface pattern 
of hnes. The screen for the optical deflection technique is 
dismantled and the pattern is illuminated by a diffusive light 
source. Due to the rectangular geometry the lines are oriented 
perpendicular to the longitudinal container axis (transverse 
mode). The deformations appearing at the arched ends die 
out after a few pattern wavelengths. The accelerometer and 
some cabling are visible on the left side of the container. The 
background of the picture is the shaker armature. 
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FIG. 


2. Sketch of the amplitude measurement technique. 



A laser beam with a diameter comparable to the wavelength 
of the surface pattern passes through a hole in the center of 
a diffusing screen and gets reflected at the oscillating fluid 
surface. In case of a line pattern the laser spreads out to a 
light streak of length s on the screen, which is recorded by 
a CCD camera. Since the light ray defining the tips of the 
streak may originate from two neighboring wave nodes the 
length s is flawed by a systematic error As of one wavelength 
27r/fc- 
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FIG. 3. Sample SI; Stationary wave amplitude 7?|max as a 
function of the reduced acceleration e = a/ac — 1 taken at the 
drive frequencies (a): Q.jl'K = 60Hz, (b): 80Hz, (c): 120Hz, 
(d) 160Hz. The maximum elevation ?7|max is derived from the 
measured slope via the relation 7)|max = 9?7/9a;|max/fcc. Solid 
symbols: experimental data, solid line: result of the weakly 
nonlinear perturbation analysis of Ref. open symbols: 

results of the full scale numerical simulation. The band be- 
tween the dotted lines in (a) indicates the systematic error 
of the measuring technique. Error bars indicate statistical 
fluctuations as obtained by a sequence of 25 individual mea- 
surements. At small e the statistical fluctuation increase since 
the length of the light streak is small. At large e the error in- 
creases again because of a beginning defect dynamics, which 
destroys the stationarity and the coherence of the line pattern. 
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driving frequency 0^271 (Hz) reduced acceleration 



FIG. 4. The coupling coefficient l/F (scaled by fc^) ob- 
tained from sample SI (S2). The squares (circles) denote 
fitted values for T as compiled from the data of Fig. ^ (Fig.^. 
The solid (dashed) line is the prediction of the perturbation 
theory (12). 



FIG. 6. Sample S2 at T = 30°C: Stationary wave slope 
9r;/9x|max at the drive frequencies Q/2tt — 80Hz (squares), 
lOOHz (circles), 140Hz (diamonds). 
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FIG. 5. Sample SI: Bifurcation diagram for the maximum 
surface slope 9??/i9a;|max. Symbols and parameters as in Fig. ^ 
The solid line is the asymptotics given by Eq. for fl ^ 00. 
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FIG. 7. Sample S'2 at dilferent temperatures: (a) Bifurca- 
tion diagram 9?7/9a;jmax at the drive frequency Q/2tt — 80Hz. 
The viscosities are f = 94cS (squares), 85cS (circles), 73cS 
(diamonds), 63cS (down triangles), 8.35cS (stars), (b) The 
coupling coefficient l/F (scaled by k^) as a function of the 
viscosity The symbols show the coefficient estimated by 
ffiting the data of (a) for e < 1%. The prediction of the 
perturbation analysis (solid line) becomes unreliable for vis- 
cosities above u ~ 50cS. 
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8.35 
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TABLE I. Fluid specifications at different temperatures for 
the low (high) viscosity sample SI (S2) 



